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FOREWORD 

This  interim  report  describes  work  performed  at  the  Columbus  Aircraft 
Division  of  Rockwell  International  during  the  period  from  1 June  1977  to 
28  February  1978.  This  work  was  supported  by  the  Office  of  Naval 
Research,  under  Contract  N00014-77-C-0271.  The  Scientific  Officer  for 
the  project  was  Dr.  R.  E.  Whitehead.  It  is  the  first  phase  of  an  effort 
to  develop  a computer  program  for  predicting  the  V/STOL  performance  of 
thrust  augmenting  ejector  aircraft. 

The  results  presented  here  owe  much  to  the  contributions  of  other 
members  of  the  technical  staff,  in  particular  J.  K.  McCullough  and 
J.  H.  DeHart.  We  would  like  to  thank  all  of  them  for  their  valuable 

suggestions. 
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ABSTRACT 


A method  has  been  developed  for  calculating  the  static  performance 
of  thrust  augmenting  ejectors  by  matching  a viscous  solution  for  the 
flow  through  the  ejector  to  an  inviscid  solution  for  the  flow  outside 
the  ejector.  In  effect,  the  ejector  shroud  is  considered  to  be  "flying" 
in  the  secondary  velocity  field  induced  by  the  entrainment  of  the 
primary  jets.  A two-dimensional  analysis  utilizing  a turbulence 
kinetic  energy  model  for  the  inner,  jet  mixing  solution  and  potential 
flow  singularities  for  the  outer,  induced  flow  is  described.  This 
approach  offers  the  advantage  of  including  external  influences  on  the 
flow  through  the  ejector.  Comparisons  with  data  are  presented  for  an 
ejector  having  a single  central  nozzle  and  Coanda  jets  on  the  walls. 

The  accuracy  of  the  matched  solution  is  found  to  be  especially  sensitive 
to  the  jet  flap  effect  of  the  flow  just  downstream  of  the  ejector  exit. 
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INTRODUCTION 


The  U.  S.  Navy  is  developing  several  categories  of  V/STOL  aircraft 
for  use  with  smaller  carriers,  as  a more  economical  means  of  maintain- 
ing sea  control.  Ejector  thrust  augmentation  offers  a solution  to  many 
of  the  problems  that  arise  in  the  design  of  such  V/STOL  aircraft.  A 
thrust  augmenting  ejector  is  a pneumatic  device  through  which  a large 
mass  of  air  drawn  from  the  atmosphere  is  pumped  by  the  entrainment  of  a 
primary  jet,  as  shown  in  Figure  1.  The  entrained  air  is  accelerated  by 
the  turbulent  mixing  of  the  two  streams  within  the  ejector.  By  Newton's 
law  of  action  and  reaction,  a force  which  is  equal  and  opposite  to  the 
momentum  change  of  the  accelerated  fluid  is  experienced  by  the  ejector. 
Thus,  significant  increases  in  the  thrust  of  turbojet  and  turbofan 
engines  can  be  obtained  by  diverting  the  exhaust  flow  through  an  ejector 


Figure  1.  Entrainment  by  the  Primary  Jet  Induces  a Secondary  Flow 
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Since  ejectors  can  be  designed  to  deflect  as  well  as  augment  the 
engine  thrust,  they  can  be  used  to  give  an  aircraft  having  a basic 
thrust  to  weight  ratio  less  than  one  the  direct  jet  lift  necessary  for 
vertical  takeoffs  and  landings.  In  effect,  the  aircraft  is  given  a 
variable  bypass  propulsion  system:  mass  flow  and  thrust  are  augmented 
by  passing  the  flow  through  an  ejector  for  vertical  flight,  and  reduced 
by  exhausting  the  flow  through  a conventional  nozzle  for  horizontal 
flight.  When  the  ejector  is  integrated  with  the  wing  to  produce  a lift- 
propulsion  system,  separate  reaction  jets  are  not  required  for  control 
during  hover;  control  forces  are  generated  by  differential  action  of  the 
ejectors.  In  addition,  the  ejector  wing  has  good  transition  and  STOL 
performance,  because  the  exhaust  flow  acts  like  a jet  flap  to  increase 
the  circulation  lift  of  the  wing. 

Although  consideration  of  the  ejector  for  aircraft  thrust  augmenta- 
tion began  more  than  thirty  years  ago,(^  development  has  been  largely 
by  experiment.  However,  analytic  methods  are  necessary  for  conceptual 
studies  and  to  reduce  test  requirements.  In  order  to  calculate  ejector 
performance  without  solving  the  full  Navier-Stokes  equations,  some 
approximations  must  be  made.  The  analyses  that  have  been  developed  are 
broadly  based  on  von  Karman's  now  classical  approach, which  utilizes 
streamwise  integration  of  the  governing  equations.  If  the  ejector  is 
relatively  long  and  the  diffuser  angle  is  small,  gradients  of  the  normal 
stress  and  the  variation  of  pressure  across  the  flow  can  be  neglected. 

This  reduces  the  governing  elliptic  equations  to  a parabolic  set  which 
can  be  solved  by  marching  through  the  ejector  in  the  streamwise  direc- 
tion. The  solution  is  obtained  by  iterating  on  the  inlet  velocity  until 
the  exhaust  pressure  matches  the  atmospheric  pressure  outside  the  ejector. 
A method  incorporating  a mixing  length  model  for  the  turbulence  was 
developed  by  Gilbert  and  Hill^)  for  a simple  ejector  having  a single 
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primary  jet,  while  DeJoode  and  Patankar^^  employed  a two-equation 
model  for  the  turbulence  in  an  ejector  with  hypermixing^ ) and  Coanda 
'lets. 

Even  though  the  elliptic  boundary  value  problem  can  thus  be  trans- 
formed to  an  initial  value  problem  which  is  more  easily  solved,  the 
basic  elliptic  character  of  the  flow  field  is  unchanged.  This  means 
that  there  are  cases  for  which  the  classical  momentum  theories  predict 
the  wrong  performance,  or  do  not  yield  a solution,  at  all.  For  instance, 
as  the  walls  of  the  ejector  shroud  are  removed  to  infinity,  the  predicted 
augmentation  ratio  does  not  reduce  to  unity,  as  it  must  in  the  limit  of 
an  isolated  turbulent  jet.  Also,  when  the  ejector  is  short  or  the  dif- 
fuser angle  is  large,  the  exhaust  pressure  is  less  than  atmospheric 
pressure.  This  pressure  difference  is  supported  by  the  momentum  of  the 
exhaust  jet,  while  the  jet  momentum  depends  in  turn  on  the  pressure  dif- 
ference. Thus,  the  exhaust  pressure  becomes  a floating  boundary  condi- 
tion, and  a unique  solution  to  the  initial  value  problem  cannot  be 
obtained  by  classical  methods. 

The  purpose  of  this  paper  is  to  present  an  ejector  analysis  not  sub- 
ject to  these  limitations.  The  primary  elliptic  effects  are  included 
by  iterating  between  a parabolic  solution  for  the  flow  through  the  ejec- 
tor and  an  elliptic  solution  for  the  flow  outside  the  ejector.  This 
technique  is  similar  to  that  used  in  coupling  a solution  for  the  dis- 
placement thickness  of  a wing  boundary  layer  to  a solution  for  the 
external  flow.  A parabolic  solution  for  the  turbulent  mixing  within 
the  ejector  was  developed  according  to  classical  momentum  theories. 

For  the  elliptical  flow  outside  the  ejector,  a solution  based  on 
Bcvilaqua's  lifting  surface  theory^)  was  developed.  Briefly,  the  aug- 
menting force  is  related  to  the  lift  on  a wing,  and  the  circulation  is 
determined  as  a function  of  the  jet  entrainment. 
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A two-dimensional  analysis  utilizing  a turbulence  kinetic  energy 
model  for  the  inner,  viscid  solution  and  vortex  lattice  methods  for 
the  outer  inviscid  solution  is  described.  Results  are  presented  for 
an  ejector  having  a single  center  nozzle  and  jets  on  the  walls, 
but  other  configurations  can  be  analyzed  with  the  same  approach.  In 
the  next  section  an  outline  of  ejector  theory  is  presented  to  introduce 
the  mathematical  models  used  in  this  analysis.  The  solution  algorithms 
and  the  method  of  iteration  are  described  in  the  following  sections. 

The  predictions  of  this  new  model  are  compared  with  classical  solutions 
and  experimental  data  in  the  final  section. 


PRINCIPLE  OF  THRUST  AUGMENTATION 


Circulation  Thrust 

Although  ejector  thrust  augmentation  may  seem  to  utilize  a new  princi- 
ple of  lift  generation,  it  actually  involves  no  more  than  a novel  appli- 
cation of  the  familiar  circulation  theorem  of  aerodynamic  lift.  An 
isolated  jet  induces  an  essentially  lateral  flow  of  entrained  air,  as 
sketched  in  Figure  2.  However,  the  distributions  of  pressure  and  vel- 
ocity in  the  flow  entrained  by  the  jet  are  altered  by  the  presence  of  a 
shroud.  A circulation  which  redirects  the  entrained  flow  through  the  ejec- 
tor is  generated  around  each  of  the  shroud  sections,  as  shown  in  Figure  3. 
The  shroud  can  therefore  be  considered  to  be  "flying"  in  the  velocity 
field  of  the  flow  entrained  by  the  jet,  and  it  experiences  a force 
analogous  to  the  lift  developed  on  a wing  fixed  in  a moving  stream. 


Figure  2.  Streamlines  of  the  Flow  Induced  by  a Free  Jet 


5 


Figure  3.  A Circulation  Redirects  the  Flow  through  the  Ejector 

According  to  this  lifting  surface  theory, the  thrust  augmentation  <£ 
can  be  defined  as  the  ratio  of  the  primary  jet  thrust  T plus  the  "lift" 
on  the  shroud  F to  the  isentropic  thrust  of  the  primary  mass  Tj: 

T + F 


Tl 


(1) 


The  thrust  augmentation  results  from  the  fact  that  the  interaction 
between  the  flow  induced  by  the  entrainment  of  the  jet  and  the  vorticity 
bound  in  the  sections  of  the  shroud  generates  a pair  of  equal  and 
opposite  forces.  The  origin  of  these  forces  can  be  understood  by  a 
consideration  of  the  interaction  between  a sink  of  strength  Q,  which 
represents  a section  of  the  jet,  and  a vortex  of  strength  T,  which 
represents  a segment  of  the  vortex  sheet  in  the  shroud.  These  singu- 
larities are  a distance  r apart,  as  shown  in  Figure  4.  At  the  vortex, 
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Figure  4.  Mutually  Induced  Forces  on  a Sink-Vortex  Pair 

the  sink  induces  a velocity  of  magnitude  Q/2^r,  directed  along  r.  The 
vortex  therefore  experiences  a force  prQ/2wr,  perpendicular  to  r.  At 
the  sink,  the  vortex  induces  a velocity  of  magnitude  r/2-rrr,  perpendicu- 
lar to  r.  The  sink  therefore  experiences  a force  pQr/27rr,  also  perpen- 
dicular to  r,  but  opposite  to  the  force  on  the  vortex.  The  net  effect 
of  the  interactions  between  all  the  sinks  and  vortices  is  a force  which 
increases  the  thrust  of  the  jet,  and  an  equal  but  opposite  reaction  on 
the  shroud. 

The  force  on  the  shroud  can  be  recognized  as  the  vortex  force  given 
by  the  Kut ta- Joukowski  theorem  for  airfoil  lift.  This  force  appears  in 
the  pressure  distribution  on  the  surface  of  the  shroud,  primarily  as  a 
leading  edge  suction.  However,  the  thrust  on  the  jet  sinks  may  not  be 
as  familiar.  This  force  is  conceptually  similar  to  the  ram  drag  that 
develops  on  an  aircraft  engine  inlet,  but  it  must  be  remembered  that  the 
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sink/vortex  interaction,  as  described,  is  only  valid  in  irrotational 
flows.  The  fluid  entrained  into  the  jet  becomes  rotational,  so 
that  an  understanding  of  how  this  force  actually  developes  requires 
consideration  of  the  jet  mixing  process. 

Jet  Thrust  Augmentation 


The  mixing  process  is  basically  an  inelastic  collision  between  the 
jet  and  the  surrounding  fluid.  As  such,  jet  mixing  is  governed  by  the 
same  laws  of  momentum  and  energy  conservation  as  simple  collisions 
between  discrete  particles.  In  free  jet  mixing,  the  momentum  flux  is 
conserved  and  there  is  a corresponding  loss  of  kinetic  energy  to  turbu- 
lence and  heat.  However,  the  circulation  accelerates  the  entrained  flow 
entering  the  ejector  and,  according  to  Bernoulli's  equation,  the  pres- 
sure drops.  Inside  the  ejector,  before  mixing,  the  velocities  in  the 
primary  and  secondary  flows  are  given  approximately  by 


(2) 


* / 2AP  n 

us  = (~r) 


1/2 


(3) 


in  which  AP  is  the  pressure  drop  induced  by  the  circulation.  Momentum 
is  conserved  during  the  mixing  process  itself,  so  that 

mp  Up*  + ms  Ug*  = (nip  + ms)  Um  (4) 

The  mixed  stream  decelerates  as  it  leaves  the  low  pressure  region,  and 
its  final  velocity  is 

/ 2AP 

um  = (Um  (5) 

The  ratio  of  the  final  momentum,  (mp  + ms)Um,  to  the  initial  momentum, 
mpUp,  may  be  evaluated  by  substituting  in  turn  for  Um,  then  Um,  and 
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finally  for  UpV  and  Ug.  Performing  the  substitutions  yields  for  this 


<t>  = Vl  + 2M  (H2  + H)  1 - H (6) 

in  which  M is  the  mass  flux  ratio,  M = mg/mp,  and  H is  the  ratio  of  the 
static  pressure  drop  to  the  initial  kinetic  energy,  H = AP/V2pUp2. 

If  the  pressure  drop  is  negligibly  small,  H<<:1  and  the  solution 
reduces  to  the  case  of  free  jet  mixing.  Momentum  is  conserved,  <j>  = 1. 

On  the  other  hand  if  the  pressure  drop  is  relatively  large,  H»1  and 
the  momentum  ratio  simplifies  to 

<t>  = (1  + M)V2  (7) 

so  that  the  thrust  increases  with  the  entrainment  ratio.  The  dependence 
of  <f>  on  intermediate  values  of  H is  sketched  in  Figure  5. 


Figure  5.  .Jet  Thrust  Augmentation  as  a Function  of  Entrainment  and 
Pressure  Drop 
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Thus,  ejector  thrust  augmentation  is  generated  by  the  turbulent  mix- 
ing of  the  primary  jet  and  secondary  stream  in  a region  of  reduced 
pressure.  The  increase  in  thrust  that  results  when  the  mixed  flow 
returns  to  ambient  pressure  appears  as  a reaction  force  on  the  ejector 
shroud.  Although  these  forces  are  related  to  the  mutually  induced 
forces  on  a vortex  and  sink  in  irrotational  flow,  the  flow  through  the 
ejector  actually  includes  regions  of  interacting  irrotational  and  turbu- 
lent fluid,  subject  to  lateral  straining,  and  streamwise  curvature,  with 
variations  of  temperature  and  density.  In  the  following  section  a 
method  of  calculating  these  forces  in  a real  fluid  will  be  developed 
from  the  principles  outlined  in  this  section. 
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INVISCID,  OUTER  SOLUTION 


Jet  Sink  Distribution 

The  Navier-Stokes  equations  for  steady  fluid  motion  are  elliptic, 
which  means  that  the  domain  of  influence  of  a point  disturbance  is  the 
entire  flow  volume.  That  is,  pressure  and  stress  gradients  transmit 
the  effect  of  local  disturbances  to  every  other  point  in  the  flow;  thus, 
the  flow  through  the  ejector  depends  on  boundary  conditions  outside  the 
ejector  as  well  as  within  the  shroud.  The  basis  of  the  present  analysis 
is  the  separation  of  the  flow  field  into  two  regions:  the  region  of 
viscous  mixing  within  the  ejector  and  the  inviscid  region  around  it. 
Significant  eLliptic  effects  are  transmitted  through  the  pressure  field 
of  the  external  flow,  while  the  turbulent  mixing  is  accurately  calcu- 
lated for  the  internal  flow.  Short  of  solving  the  full  Navier-Stokes 
equations,  matching  the  solutions  in  these  regions  is  the  only  way  to 
predict  ejector  performance.  Since  it  reduces  the  requirements  for 
computer  time  and  storage,  the  matching  procedure  has  some  advantage 
over  the  full  solution. 

The  circulation  generated  around  each  section  of  the  ejector  shroud 
will  be  calculated  by  solving  a system  of  equations  which  specify  that 
the  shroud  must  be  a streamline  of  the  flow  induced  by  the  entrainment 
of  the  jets.  For  this  calculation,  the  strengths  of  the  sinks  which 
represent  each  jet  are  considered  to  be  known.  These  strengths  are 
determined  from  the  turbulent  mixing  of  the  jets  computed  in  the  vis- 
cous solution  during  the  previous  iteration. 

The  entrainment  of  the  central  jet  is  represented  by  a series  of 
overlapping,  triangular  sink  distributions  on  the  axis  of  the  ejector, 
as  shown  in  Figure  6.  Each  distribution  is  identified  by  the  index  of 
its  central  point,  where  the  strength  of  the  sink  is  qj;  every  panel  is 
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Figure  6.  Triangular  Sink  Panels  on  the  Axis  of  the  Ejector 


the  same  length,  2s.  Tn  incompressible  flow,  the  horizontal  and  verti- 
cal components  of  velocity  induced  at  an  arbitrary  point  P(x,  y)  by 
such  a distribution  are  given  by: 

u(x,y)  = q j ^ ^2  (1  1 d-  (8) 

j 

v(x,y)  = qj  J -^2  (1  ± I*  d-  (9) 

j 

in  which  i is  a coordinate  on  the  path  of  integration  over  the  panel, 
and  r is  the  distance  [(x  - £)2  + y2]  • The  plus  sign  (+)  is  used  for 

the  rising  side  of  the  triangle,  from  -s  to  the  center,  and  the  minus 
sign  (-)  for  the  falling  side,  from  the  center  to  +s.  The  integration 
of  these  expressions  is  given  in  Appendix  I. 
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The  surface  of  the  ejector  shroud  is  represented  by  m source  panels 
of  different  lengths,  ij , and  uniform  strength,  q j , as  shown  in  Figure 
7.  Because  a single  sheet  of  sinks  cannot  provide  the  jump  in  entrain- 
ment necessary  to  model  the  presence  of  a wall  jet  on  the  inner  surface 
of  the  shroud  only,  both  the  inner  and  outer  surfaces  must  be  represented 
by  source/sink  panels.  The  horizontal  and  vertical  components  of  veloc- 
ity induced  at  an  arbitrary  point  by  a uniform  source  distribution  are 
given  by: 


u(x,y) 


(10) 


v(x,y) 


(11) 


f igurt  7.  Uniform  Source/Sink  Panels  on  the  Surface  of  the  Shroud 
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These  expressions  are  also  integrated  in  Appendix  I.  The  velocity  com- 
ponents induced  at  each  point  in  the  flow  are  obtained  by  adding  up 
the  contributions  of  every  panel. 


Shroud  Vorticity  Distribution 

A vortex  lattice  method  was  used  to  determine  the  circulation  density 
of  the  shroud.  The  continuous  vorticity  distribution  is  replaced  by  n 
discrete  vortices  of  strength  Tj  , located  at  Xj  , the  quarter  chord  of 
the  panels  shown  in  Figure  8.  It  was  found  that  better  results  were 
obtained  if  the  vortex  sheet  is  placed  on  the  inner  surface  of  the 
shroud,  rather  than  on  the  mean  camber  line.  All  the  flow  singularities 
which  represent  the  shroud  geometry  and  jet  effects  must  induce  the 
known  entrainment  (inflow)  velocities  on  the  surface  of  the  shroud. 
However,  the  source/sink  distribution  on  each  surface  already  satisfies 
this  boundary  condition.  Therefore,  the  resultant  of  the  velocities 
induced  by  all  the  other  singularities  must  be  tangent  to  the  inner 


Figure  8.  Discrete  Vortex  Lattice  on  the  Inner  Surface  of  the  Shroud 
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surface  of  the  shroud;  that  is,  the  normal  velocity  induced  by  the 
vortex  sheet  must  be  equal  but  opposite  to  the  normal  velocity  induced 
by  the  central  jet  and  opposite  wall  jet. 


The  vortex  strengths  are  determined  by  satisfying  this  boundary 
condition  at  n points  which  correspond  to  the  three  quarter  chord 
station  on  each  panel.  In  calculating  the  velocity  induced  by  the 
vortex  sheets,  it  is  convenient  to  simultaneously  consider  the  influence 
of  each  vortex  and  its  image  on  the  opposite  side  of  the  ejector.  The 
horizontal  and  vertical  components  of  velocity  induced  at  a point 
p(xi*  yj.)  by  the  vortex  pair  of  strength  I'j  at  the  points  P(xj,  y j > and 
P(xj , -yj)  are 


u(xif  yt)  = 


nzJOL  tun  )r 

2 TTTl  . 2 _ 4 2 j j 

1J  1 » J ' 


v(xi,  yi)  = 2 + ^ 2 ri 

\ 27rri j 2 2«it_ j2/  J 


2 2-i /z 

in  which  r-j  = [(x^  - Xj ) + (y^  — y j ) J is  the  distance  between 

points.  Thus,  the  contribution  to  the  velocity  normal  to  panel-i  by 


both  vortex  sheets  is 


Vv  = zC  <bij  cos  ai  * aij  sin  ai>  Pj 


in  which  the  influence  coefficients,  ajj  and  b^j , have  the  form  given 
in  Equations  (12)  and  (13),  and  is  the  angle  of  panel-i  relative  to 
the  ejector  axis.  Similarly,  the  normal  velocity  induced  by  the  central 
jet  and  the  opposite  wall  jet  is 

vs  = £ <dij  cos  ocj  - c.j  sin  04)  q (15) 
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in  which  the  range  of  the  index  j is  over  both  jets,  so  that  the 
influence  coefficients  are  given  by  Equations  (8)  and  (9)  or  Equations 
(10)  and  (11),  as  appropriate.  As  previously  noted,  the  sink  strengths 
are  known,  having  been  determined  by  the  solution  of  the  turbulent  mix- 
ing problem  in  the  previous  iteration. 

Since  the  resultant  of  the  normal  velocities  induced  by  the  jet  and 
vortex  sheet  must  be  zero,  Vv  is  set  equal  but  opposite  to  V s at  each 
control  point, 


-V 


si 


A«j  q 


(16) 


The  summation  convention  for  repeated  indices  is  intended  to  apply. 

Thus,  this  expression  represents  a set  of  simultaneous  algebraic  equa- 
tions for  the  I\.  The  influence  coefficients,  A^j,  have  the  form  given 
in  Equation  (14).  This  equation  was  solved  for  the  I'j  by  triangulariza- 
tion  of  the  coefficient  matrix. 


I 

I 


VISCOUS,  INNER  SOLUTION 

Governing  Equations 

The  entrainment  which  represents  the  effect  of  the  jets  is  calculated 
f rom  a solution  for  the  turbulent  mixing  within  the  ejector.  It  is 
possible  to  calculate  the  rate  of  entrainment  without  solving  the  com- 
plete three-dimensional  mixing  problem,  by  taking  advantage  of  the  flow 
geometry.  Since  there  is  a primary  direction  of  flow  (through  the 
ejector)  it  is  assumed  that  the  thin  shear  layer  approximation  can  be 
applied.  This  reduces  the  governing  elliptic  equations  to  a parabolic 
set.  The  effect  of  a disturbance  is  confined  to  regions  downstream  of 
the  disturbance,  and  so  the  equations  can  be  solved  by  marching  through 
the  ejector  in  the  streamwise  direction.  This  results  in  considerable 
savings  in  computer  storage  and  time  compared  to  a solution  of  the 
elliptic  mixing  problem. 


The  two-dimensional  flow  in  an  ejector  cross  section  is  governed  by 
the  continuity  equation  and  the  streamwise  momentum  equation.  The  thin 
shear  layer  approximation  means  that  the  gradients  of  the  normal  stress 
are  negligible,  and  the  pressure  P is  constant  in  each  plane  normal  to 
the  direction  of  flow.  Thus,  only  shear  stresses  caused  by  velocity 
gradients  across  the  flow  are  significant.  An  additional  assumption 
that  the  fluid  density  p is  uniform  was  also  made.  Under  tnese  assump- 
tions, the  equation  for  the  conservation  of  mass  and  momentum  through 
the  ejector  become: 


£-» 

(17) 

= dr  _ dP 

3y  dx 

(18) 
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Here,  u is  the  time  averaged  velocity  in  the  streamwise  direction,  and 
T is  the  turbulent  shear  stress.  Laminar  stresses  are  assumed  to  be 
negligible. 

In  order  to  provide  accurate  calculations  of  the  turbulent  stresses 
in  each  region  of  the  flow  (initial  and  developed  sections  of  the  free 
jet,  inner  and  outer  layers  of  the  wall  jet,  and  the  merged  region)  the 
advanced  two  equation  turbulence  model  of  Launder  and  Spalding^^  was 
used  for  turbulence  closure.  This  model  employs  differential  equations 
for  the  turbulence  kinetic  energy  and  its  rate  of  dissipation;  the  eddy 
viscosity  at  each  point  is  calculated  from  these  quantities.  This 
approach  has  previously  been  used  in  more  complex  jet  flows  by  DeJoode 
and  Patankar,^)  where  it  yielded  satisfactory  results. 


According  to  the  usual  eddy  viscosity  assumption,  the  turbulent 
stress  is  first  expressed  in  terms  of  a turbulent  viscosity  and  the 
velocity  gradient  in  the  cross  stream  direction: 


(19) 


The  two-equation  turbulence  model  developed  by  Launder  and  Spalding^) 
gives  the  turbulent  viscosity  in  terms  of  two  parameters,  for  which  two 
differential  equations  are  solved.  The  expression  for  turbulent  vis- 
cosity is: 


Cp  p k2 

( 


(20) 


where  c^  is  a constant,  k is  the  kinetic  energy  of  turbulence,  and  € 
is  the  rate  of  its  dissipation.  In  two-dimensional  parabolic  flows, 
the  governing  equations  for  k and  f are: 


Pu 


PU 


0€ 

r!x 


3k  : 
3x 

. _3_ 
dy 


_3_ 

3y 

(^#)+G  - p€ 

(21) 

ae 

f)+  <-iG  - o2  pO  f 

(22) 
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The  quantity  G is  the  rate  of  generation  of  k by  the  action  of  the 
velocity  gradient  and  is  expressed  as: 


G = 


(23) 


This  turbulence  model  involves  five  empirical  constants.  According 
to  the  recommendation  of  Launder  and  Spalding^®^  the  following  values 
of  the  constants  are  appropriate  for  ordinary  jets: 

Cg  C]^  C2  <7^  Of. 

0.09  1.44  1.92  1.0  1.3 

The  effect  of  hypermixing^ ^ can  be  simulated  by  increasing  the  value 
of  Cp  by  a factor  of  5 to  6 in  the  central  jet. 


Method  of  Solution 


The  procedure  used  to  solve  the  governing  equations  is  very  similar 

(a) 

to  the  method  devised  by  Patankar  and  Spalding.  ' The  equations  are 
transformed  from  a Cartesian  coordinate  system  (x,y)  to  a (x,o>)  system 
where  u>  is  a dimensionless  stream  function  defined  as: 

w = “t — — (24) 

^total 

where  dij,  = pudy  and  is  the  total  mass  flow  rate  in  the  computa- 

tional domain.  The  method  employs  a finite-difference  marching  proced- 
ure; from  known  conditions  at  an  upstream  cross  section,  x,  the  flow 
field  at  the  downstream  cross  section,  x + Ax,  is  computed.  This 
marching  process  is  continued  until  the  domain  of  interest  is  covered. 
The  initial  conditions  are  determined  from  the  velocities  induced  in 
the  ejector  inlet  by  the  vorticity  distribution  obtained  in  the  outer 
solution  on  the  previous  iteration.  The  finite-difference  equations 
are  formed  by  integrating  the  differential  equations  over  a small 
control  volume  surrounding  each  grid  point.  The  resulting  non-linear 
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equations  are  linearized  by  using  upstream  values  of  the  flow  variables 
to  evaluate  coefficients  involving  cross  stream  convection  and  diffusion. 
The  equations  are  solved  by  the  use  of  the  tri-diagonal  matrix  algorithm. 

The  procedure  employed  for  evaluation  of  the  unknown  pressure  grad- 
ient that  appears  in  the  momentum  equation  requires  special  mention.  A 
technique  developed  by  Sparrow,  et  al^®)  is  USed  to  determine  the 
correct  value  of  dP/dx.  Briefly,  a non-linear  equation  for  dP/dx  is 
derived  by  combining  the  momentum  equation  with  the  constraint  that  the 
value  of  dP/dx  must  result  in  values  of  u which  satisfy  the  requirements 
of  continuity.  This  equation  is  solved  by  an  iteration  algorithm  that 
is  an  adaptation  of  the  Newton-Raphson  procedure.  The  convergence  is 
fast  and  the  procedure  eliminates  the  arbitrariness  and  the  resultant 
inaccuracy  of  the  technique  originally  developed  by  Patankar  and 
Spalding. 

Jet  Entrainment  Rate 

The  jet  entrainment  is  derived  from  the  solution  for  the  turbulent 
mixing  within  the  ejector.  In  principle,  the  entrainment  can  be 
obtained  by  integrating  the  jet  velocity  profiles,  and  finding  the 
difference  in  mass  flux  between  successive  stations.  However,  practical 
difficulties  in  defining  the  jet  boundaries  introduce  considerable  error 
in  calculating  the  mass  flux,  and  this  procedure  cannot  be  used  at  all 
beyond  the  point  where  the  jets  merge.  For  these  reasons  a method  was 
developed  for  determining  the  entrainment  indirectly,  from  the  pressure 
rise  within  the  ejector. 

The  connection  between  the  entrainment  rate  and  the  pressure  rise  is 
illustrated  in  Figure  9,  which  shows  the  spreading  of  a single  jet  in  a 
straight  duct.  As  the  jet  spreads,  fluid  is  entrained  into  the  primary 
stream  and  the  velocity  of  the  secondary  stream  is  correspondingly 
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Figure  9.  Spreading  of  a Turbulent  Jet  in  a Duct 


reduced.  According  to  Bernoulli's  equation,  the  static  pressure  in  the 
secondary  stream  rises.  The  boundary  of  the  jet  continuously  adjusts 
to  equalize  the  pressure  across  the  duct  while  providing  the  necessary 
entrainment.  Mixing  continues  downstream  of  the  point  where  the  jet 
has  spread  to  the  walls,  and  the  pressure  rises  accordingly.  The 
static  pressure  can  be  related  to  an  equivalent  secondary  velocity  in 
this  region.  This  velocity  is  shown  in  Figure  9. 

The  entrainment  of  the  jet  between  two  successive  stations  can  there- 
fore be  calculated  as  the  quantity  of  fluid  that  must  be  removed  from 
the  secondary  stream  to  produce  the  observed  pressure  rise: 

me  = p(2APi.1/P)V2  A._ L - p(2APi/p)V2Ai  (25) 


in  which  Bernoulli's  equation  is  used  to  express  the  secondary  velocity 
in  terms  of  the  local  static  pressure  difference,  AP^  = Pa  - P^.  For 
the  general  case,  the  cross  sectional  area  of  the  duct,  A^,  may  vary 
from  station  to  station.  If  there  is  more  than  one  jet  within  the  duct, 
as  in  the  present  configuration,  the  total  entrainment  is  apportioned 
between  the  jets  according  to  the  relative  spreading  rates. 

The  sink  strengths  that  will  represent  the  effect  of  the  jets  in 
the  inviscid  calculation  are  determined  from  the  entrainment  of  each 
jet.  An  entrainment  velocity,  Vef  is  derived  from  the  mass  entrained 
between  successive  stations,  me,  according  to  the  definition, 

ve  = me/  p Ax  (26  ) 

in  which  Ax  is  the  distance  between  stations.  For  the  central  jet,  the 
strengths  of  the  qj  are  determined  by  setting  the  velocity  induced  at 
the  midpoint  of  each  triangular  distribution  equal  to  the  entrainment 
velocity  at  that  point.  It  is  shown  in  Appendix  II  that  qj  = -2Ve  at 
each  control  point. 

The  strengths  of  the  qj  for  the  wall  jets  on  the  surface  of  the  ejector 
shroud  are  determined  by  simultaneously  satisfying  the  known  entrainment 
(inflow)  boundary  condition  due  to  the  wall  jet  on  the  inner  surface  of  the 
shroud,  and  the  condition  of  zero  flow  through  the  outer  surface.  By 
a transformation  of  coordinates,  Equations  (10)  and  (11)  are  used  to 
calculate  the  normal  component  of  velocity,  V,  induced  by  the  j-th 
source  distribution  at  the  midpoint  of  the  i-th  panel.  In  Appendix  II, 
the  velocity  induced  at  the  midpoint  of  panel-i  by  the  source  distribu- 
tion on  panel-i  is  shown  to  be  = q^/2.  The  contribution  of  every 
panel,  including  panel-i,  is  added  at  each  control  point.  The  normal 
velocity  component  at  each  control  point  is  given  by  the  summation, 

Vi  = S (vi j cos  #ij  " uij  sin  tfij) 

j 

in  which  is  the  angle  between  the  normals  to  panels  i and  j. 
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The  induced  velocity  is  set  equal  to  the  required  entrainment  velocity 
Vc  at  each  control  point,  yielding  a set  of  m simultaneous  algebraic 
equations  for  the  q^ , 


in  which  the  influence  coefficients  give  the  velocity  induced  at  the 
point  i by  a source  of  unit  strength  on  panel  j.  This  matrix  equation 
is  solved  for  the  qj  by  triangularizat ion  of  the  coefficient  matrix, 

The  solution  yields  sinks  (negative  sources)  on  the  inner  surface 
of  the  shroud  and  positive  sources  on  the  outer  surface  of  the  shroud. 

Sine  the  reaction  to  the  thrust  of  the  mixed  stream  appears  in  the 
pressure  distribution  on  the  surface  of  the  shroud,  the  whole  effect  of 
the  turbulent  mixing  is  represented  by  the  strength  of  these  sinks. 

The  velocities  induced  by  the  sinks  determine  the  strength  of  the  vortex 
sheet  and  the  pressures  on  the  shroud.  Conversely,  the  requirement  that 
the  total  vorticity  must  satisfy  the  Kutta  condition  determines  the  rate 
of  entrainment  by  controlling  the  initial  velocity  of  the  secondary 
stream.  The  actual  strength  of  the  sinks  and  the  magnitude  of  the  cir- 
culation are  determined  by  iteration  between  the  turbulent  and  inviscid 
solutions,  as  described  in  the  next  section. 
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SOLUTION  MATCHING  PROCEDURE 

Method  of  Iteration 

The  inner  and  outer  routines  are  incorporated  into  a single  computer 
program  which  yields  a matched  solution  by  iterating  between  the  rate 
of  entrainment  by  the  jets  and  the  circulation  around  the  shroud.  For 
a simple  ejector  having  a single  primary  jet,  the  method  of  iteration 
is  relatively  straightforward.  The  sequence  is  outlined  in  the  following: 

1.  An  assumption  is  made  for  the  magnitude  of  the  secondary  velocity 
at  the  inlet  of  the  ejector. 

2.  The  entrainment  of  the  jet  is  computed  from  this  initial  condi- 
tion by  marching  through  the  ejector  with  the  turbulent  mixing  routine. 

3.  The  equivalent  sink  strengths  within  the  ejector  are  calculated 
from  the  pressure  gradient.  Downstream  of  the  ejector  exit,  it  is 
specified  that  the  sink  strengths  decay  as  in  a free  jet;  that  is,  as 


4.  The  vorticity  distribution  is  computed  with  the  inviscid  routine, 
by  requiring  that  the  shroud  be  a streamline  of  the  flow  induced  by  the 
sink  distribution. 

5.  The  velocity  induced  at  the  ejector  inlet  by  all  the  sinks  and 
vortices  is  calculated  and  compared  to  the  initial  guess. 

If  these  agree,  the  solutions  are  matched.  Otherwise,  steps  2 through 
5 are  repeated  until  convergence  is  achieved,  or  the  routines  break  down 
due,  for  example,  to  flow  separation  from  the  inner  surfaces  of  the  shroud. 

Since  such  ejector  configurations  separate  at  typically  low  inlet  and 
diffuser  area  ratios,  without  producing  significant  augmentation,  these 
are  of  limited  interest.  When  wall  jets  are  added  to  prevent  separation, 
the  maximum  augmentation  increases,  but  a new  flow  phenomenom  complicates 
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the  calculation  procedure.  Curvature  of  the  jet  sheet  leaving  the 
trailing  edge  of  the  flap,  as  shown  in  Figure  10,  supports  a Low  pres- 
ure  region  behind  the  ejector.  Morel  and  Lissaman^^  pointed  out  the 
similarity  of  this  action  to  that  of  the  jet  flap,  and  described  the 
phenomenon  as  a "jet  flap  diffuser." 

Because  the  pressure  difference  is  supported  by  the  momentum  of  the 
exhaust  jet,  while  the  jet  momentum  depends  in  turn  on  the  pressures  at 
the  ejector  exit,  another  calculation  is  required  to  determine  the 
influence  of  the  jet  flap  diffuser.  This  influence  was  determined  using 
an  approach  developed  from  the  classical  jet  flap  theory  of  Spence. 

The  pressure  difference  across  the  trailing  jet  sheet  is  related  to  the 
strength  of  an  equivalent  vortex  sheet, 

puy  = AP  (29) 


T 


so  that  the  basic  mathematical  problem  becomes  finding  a vorticity  dis- 
tribution which  makes  the  jet  sheet  a streamline  of  the  flow.  Since  the 
pressure  difference  is  balanced  by  inertia  forces  due  to  curvature  of 
the  jet  sheet,  as  shown  in  Figure  10,  the  radius  of  jet  curvature,  R, 
is  given  by 

| = AP  (30) 

in  which  J is  the  thrust  of  the  wall  jet  at  the  ejector  exit.  To  a 
first  approximation,  both  the  jet  thrust  and  radius  of  curvature  can 
be  assumed  constant. 

These  two  additional  boundary  conditions  for  the  shape  and  strength 
of  the  jet  flap  diffuser  are  satisfied  as  part  of  the  iteration  to  match 
the  inner  and  outer  solutions.  The  calculation  sequence  in  this  case  is 
similar  to  that  outlined  previously: 

1.  An  assumption  is  made  for  the  secondary  velocity  at  the  ejector 
inlet. 

2.  The  turbulent  mixing  routine  is  used  to  compute  the  development 
of  the  internal  flow  up  to  the  ejector  exit. 

3.  The  curvature  of  the  jet  sheet  which  defines  the  boundary  of  the 
jet  flap  diffuser  is  calculated  according  to  Equation  (30)  from  the  exit 
pressure  and  wall  jet  thrust. 

4.  Considering  this  jet  sheet  to  represent  a solid  boundary,  compu- 
tation of  the  turbulent  mixing  is  continued  through  the  jet  flap  diffuser. 

5.  The  equivalent  sink  strengths  are  determined  within  the  ejector 
and  jet  flap  diffuser. 

6.  The  vorticity  distribution  on  both  the  shroud  and  jet  flap  dif- 
fuser are  computed  with  the  inviscid  routine. 

7.  The  velocity  induced  at  the  ejector  inlet  by  all  the  sinks  and 
vortices  is  calculated  and  compared  to  the  initial  guess.  Steps  2 
through  7 are  repeated  until  these  agree. 
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In  this  calculation  the  length  of  the  jet  flap  diffuser  is  defined 
to  he  the  point  where  the  jet  sheets  become  parallel  to  each  other  and 
the  axis  of  the  ejector.  When  the  iteration  converges  and  the  solutions 
are  matched,  the  pressure  within  the  ejector  reaches  atmospheric  pres- 
sure at  this  point.  In  effect,  the  Kutta  condition  for  the  vorticity 
distribution  is  satisfied  at  the  end  of  the  jet  flap  diffuser,  rather 
than  at  the  trailing  edge  of  the  ejector  shroud.  In  general,  this 

approximation  gives  good  results.  However,  if  the  jets  do  not  merge 

(13) 

until  far  downstream  of  the  ejector  exit,  Marstersv  ' has  observed 
that  the  jets  can  curve  back  toward  the  axis  and  the  pressures  within 
the  jet  flap  diffuser  can  actually  rise  above  atmospheric  pressure. 
Further  development  of  the  method  will  be  required  in  order  to  obtain 
a solution  in  such  cases. 

Evaluation  of  the  Thrust  Augmentation 

The  thrust  of  the  ejector  is  equal  to  the  sum  of  the  initial  jet 
thrust  plus  the  force  on  the  shroud.  In  principle,  it  can  be  evaluated 
either  by  integrating  the  thrust  of  the  mixed  flow  at  the  ejector  exit, 
or  by  finding  the  resultant  of  the  pressure  and  shear  stresses  on  the 
surface  of  the  shroud.  In  practice,  integrating  the  stream  thrust  is 
relatively  easy  to  do,  while  a large  number  of  panels  and  an  additional 
calculation  of  the  skin  friction  are  required  to  evaluate  the  force  on 
shroud.  However,  because  the  pressure  distribution  is  needed  to  estimate 
surface  loads  and  hinge  moments,  both  calculations  are  performed. 


The  stream  thrust  is  evaluated  at  the  exit  of  the  ejector.  Although 
any  station  in  the  jet  wake  should  give  the  same  result,  the  somewhat 
artificial  treatment  of  the  flow  in  the  jet  flap  diffuser  led  to  selec- 
tion of  the  exit  plane  for  this  calculation.  The  thrust  of  the  ejector 
is  given  by 


r 


(Pa  - Pe)Ae 


(31) 
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in  which  Pe  and  Ae  are  the  static  pressure  and  area  at  the  exit. 
Because  the  flow  velocities  are  constant  over  a small  control  volume 
surrounding  each  grid  point,  integration  of  the  stream  thrust  involves 
a simple  summation  of  the  thrust  increment  from  each  control  volume. 
The  static  pressure  is  constant  across  the  exit.  It  should  be  noted 
that  even  though  the  pressure  force  is  negative,  lowering  the  exhaust 
pressure,  as  with  the  jet  flap  diffuser,  results  in  a net  thrust 
increase.  This  is  because  the  momentum  flux  is  increased  more  than 
the  pressure  force  is  reduced. 


The  surface  pressure  distribution  is  determined  from  the  tangential 
velocities  induced  at  the  surface  of  each  panel.  Actually,  these  pres- 
sures are  more  representative  of  the  pressures  at  the  outer  edge  of 
the  wall  jets.  The  high  curvature  of  the  wall  jets  over  the  leading 
edge  of  the  shroud  induces  a significant  radial  pressure  gradient 
across  the  jet.  The  resulting  decrease  in  pressure  at  the  surface  is 
given  by 


in  which  puj/Rw  is  the  inertia  force  which  resists  the  turning  of  the 
jet,  and  o is  the  jet  thickness.  This  contribution  to  the  surface  pres- 
sure has  not  been  included  in  the  present  analysis. 
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RESULTS  AND  DISCUSSION 


Comparison  with  Experiment 

In  order  to  evaluate  the  basic  lifting  surface  theory  that  the  force 
on  the  shroud  is  related  to  the  lift  on  a wing,  the  prediction  of  this 
analysis  will  be  compared  to  experimental  data.  There  are  no  exact 
solutions  for  the  surface  pressures  or  thrust  augmentation;  however, 
comparisons  will  be  made  with  the  results  of  classical  analyses  to 
determine  if  the  primary  elliptic  effects  have  been  correctly  predicted. 
A sketch  of  the  test  ejector  is  shown  in  Figure  11.  It  combines  a 
single  central  nozzle  with  Coanda  jets  on  the  inner  surface  of  the 
shroud.  A slot  in  each  endwall  at  the  ejector  throat  provides  a 
boundary  layer  control  jet  to  prevent  separation  of  the  flow  from 
these  surfaces.  In  this  configuration,  607,  of  the  primary  flow  is  in 


the  central  jet,  177„  goes  to  each  of  the  Coanda  jets,  and  the  remaining 


Figure  11.  Sketch  of  the  Experimental  Ejector 
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b “ of  the  primary  flow  is  divided  between  the  two  endwall  jets.  The 
ejector  has  a span  of  36  cm,  a length  of  27  cm,  and  has  a throat  12  cm 
wide.  The  inlet  area  ratio  is  approximately  11. 


The  test  ejector  is  supported  on  four  flexible  rods  which  are 
attached  to  a supporting  frame.  The  net  thrust  is  measured  with  two 
load  cells  installed  on  connecting  arms  between  the  frame  and  ejector. 

High  pressure  air  is  supplied  to  the  primary  jets  from  a single  source, 
but  the  rate  of  flow  to  each  nozzle  is  measured  separately,  with  a 
calibrated  venturi  in  each  of  the  feed  lines.  Internal  screens  and 
baffles  are  used  to  smooth  out  spanwise  variations  in  the  je  velocity, 
so  that  the  stagnation  pressure  is  measured  with  a total  pre  .sure  probe 
in  the  exit  of  each  nozzle.  The  isentropic  reference  thrust  is  calcu- 
lated as  the  product  of  the  measured  primary  mass  flow  and  the  velocity 
defined  by  an  isentropic  expansion  from  the  measured  stagnation  pres- 
sure to  atmospheric  pressure.  It  is  estimated  that  the  error  in  these 
measurements  is  within  ±37.  of  the  mean. 

In  Figure  12  the  calculated  change  in  the  thrust  augmentation  ratio 
with  the  diffuser  area  ratio  is  compared  to  experimental  values.  At  low 
diffuser  area  ratios  the  thrust  augmentation  is  underpredicted  by  approxi- 
mately 67.,  while  good  predictions  of  the  maximum  augmentation  are 
obtained.  This  result  is  a consequence  of  the  approach  taken  in  calcu- 
lating the  jet  flap  effect.  Because  the  length  of  the  jet  flap  diffuser 
is  defined  by  the  point  where  the  jet  sheets  become  parallel,  the  jet 
diffuser  length  goes  to  zero  as  the  diffuser  angle  of  the  duct  is 
reduced.  Most  of  the  discrepancy  at  the  low  diffuser  area  ratios  can 
probably  be  attributed  to  this  effect.  With  this  perspective,  the 
agreement  between  analysis  and  experiment  can  be  judged  satisfactory. 

It  should  be  noted  that  a plane  slot  nozzle  was  installed  in  the  test 
ejector  to  provide  a simple,  two-dimensional  jet  for  these  comparisons. 

The  augmentation  can  be  significantly  increased  by  installing  a hyper- 
mixing or  multi-lobe  nozzle. 
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The  calculated  jet  boundaries  are  compared  with  the  measured 
boundaries  in  Figure  13.  Since  the  turbulence  constants  were  not 
adjusted  for  this  case,  but  derived  from  other  flows,  the  agreement 
is  particularly  good.  The  predictions  of  the  shape  and  length  of  the 
jet  flap  diffuser  are  also  satisfactory.  In  Figure  1A  the  calculated 
velocity  distributions  at  three  stations  within  the  ejector  are  com- 
pared. The  profiles  from  the  inner,  viscous  solution  show  the  spread- 
ing of  the  jets,  as  well  as  the  reduction  in  secondary  velocities. 
Because  the  diffuser  walls  are  diverging  in  this  case,  not  all  the 
change  in  secondary  velocity  is  due  to  entrainment;  both  the  jet  and 
econdary  stream  decelerate  as  a result  of  the  change  in  cross  sectional 
a rea . 
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Figure  13.  Comparison  of  Calculated  (--)  and  Measured  (•)  Jet  Spreading  Kates 


Figure  14.  Comparison  of  Velocity  Distributions  Calculated  by  Viscous 
and  Inviscid  Solutions 


Due  to  the  assumption  that  the  static  pressures  are  constant  at  each 
axial  station,  the  secondary  velocities  in  the  viscous  solution  are 
uniform;  further,  there  is  no  transverse  velocity  component.  The 
inviscid  velocity  distributions,  shown  on  the  other  side  of  the  ejector, 
indicate  the  extent  of  the  actual  skewness  and  the  magnitude  of  the 
transverse  velocities.  Since  the  jets  are  replaced  by  equivalent  sinks 
in  the  inviscid  solution,  the  jet  profiles  are  not  seen  in  this  case. 

The  secondary  velocities  are  generally  uniform  wi thin  the  diffuser 
section.  However,  the  effect  of  the  high  inlet  curvature  reduces  the 
secondary  mass  flow  through  the  ejector:  since  the  inlet  pressures  are 
matched,  the  speed  of  the  inlet  flow  is  the  same  in  each  solution,  but 
the  inclination  of  the  inviscid  velocity  vectors  produces  a "tilt  loss" 
in  axial  mass  flow.  Curvature  effects  can  be  added  to  the  viscous  solu- 
tion, and  this  will  provide  a better  match  of  the  mass  flows.  The  pre- 
dicted pressure  rise  through  the  ejector  and  the  sink  strengths  derived 
from  it  are  shown  in  Figure  15. 


Figure  15.  Predicted  Static  Pressure  Gradient  and  Equivalent  Sink  Strengths 
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Prediction  of  Elliptic  Effects 

According  to  classical  momentum  theories , ^ the  thrust  augmenta- 
tion ratio  increases  monotonically  with  ejector  inlet  area  ratio.  In 
the  limit  of  very  large  inlet  area  ratios,  the  augmentation  approaches 
<t>  = 2.0,  with  no  diffusion.  With  diffusion,  correspondingly  larger 
values  of  augmentation  are  predicted.  This  is  because  the  parabolic 
flow  assumption  constrains  the  secondary  flow  to  be  always  parallel  to 
the  ejector  axis.  The  lifting  surface  theory^)  correctly  predicts  that 
the  augmentation  ratio  approaches  unity  in  this  limit,  as  it  must  for 
an  isolated  turbulent  jet.  But,  as  the  inlet  area  ratio  approaches 
unity  in  the  opposite  limit,  an  unrealistic  increase  in  the  augmenta- 
tion is  predicted.  This  is  because  the  effect  of  the  circulation 
velocities  in  reducing  the  rate  of  entrainment  are  neglected. 


Because  it  represents  a merging  of  these  two  theories,  the  present 
interaction  analysis  has  the  correct  behavior  in  both  limits.  Figure 
16  shows  the  predicted  variation  of  the  thrust  augmentation  ratio  as  a 


Figure  16.  Predicted  Effect  of  Inlet  Area  Ratio 
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function  of  inlet  area  ratio,  for  a constant  diffuser  area  ratio  of  1.8. 
The  length  of  the  shroud  was  kept  constant,  so  that  the  ejector  becomes 
relatively  long  at  low  inlet  area  ratios.  In  this  case  the  parabolic 
flow  approximation  is  valid,  and  the  augmentation  is  seen  to  initially 
increase  with  the  inlet  area  ratio.  This  is  as  predicted  by  the 
momentum  theories.  As  the  sides  of  the  shroud  are  moved  farther  apart, 
the  strength  and  influence  of  the  circulation  is  diminished,  and  the 
augmentation  begins  to  decrease.  This  is  according  to  the  lifting 
surface  theory.  Thus,  the  correct  behavior  has  been  predicted  in  each 
l imi t . 


35 


1 


I 


CONCLUSIONS 

A viscous/inviscid  interaction  analysis  has  been  used  to  extend 
classical  momentum  theories  of  ejector  thrust  augmentation.  The 
primary  elliptic  effects  have  been  included  by  iterating  between  a 
parabolic  solution  for  the  flow  through  the  ejector  and  an  elliptic 
solution  for  the  flow  outside  the  ejector.  Briefly,  a calculation  of 
the  rate  of  entrainment  by  the  turbulent  jets  is  used  to  determine  the 
equivalent  sink  strengths.  The  requirement  that  the  ejector  shroud 

I must  be  a streamline  of  the  flow  induced  by  these  sinks  is  then  used 

to  evaluate  the  circulation  generated  around  the  shroud.  The  influence 
of  the  circulation  is  included  in  the  next  iteration  for  the  rate  of 
entrainment.  Comparison  of  the  calculated  thrust  augmentation  with 
experimental  data  establishes  confidence  in  the  ability  to  predict  the 
complex  ejector  flowfield  with  this  approach.  In  addition,  greater 
understanding  of  the  principle  of  ejector  thrust  augmentation  is  obtained 
from  the  analysis. 

- 

The  present  analysis  can  be  improved  by  including  compressibility 
and  jet  temperature  effects,  and  by  refining  the  representation  of  the 
flow  singularities.  More  importantly,  curvature  effects  should  be 
added  to  both  the  turbulence  equations  and  the  mean  flow  equations  in 

I 

the  viscous  solution,  and  additional  development  of  the  jet  flap  dif- 
fuser scheme  should  be  undertaken. 
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APPENDIX  I 

ELEMENTARY  SINK  DISTRIBUTIONS 


In  this  Appendix  the  methods  of  potential  flow  theory  are  used  to 
derive  expressions  for  the  velocity  induced  by  elementary  sink  dis- 
tributions. This  derivation  is  based  on  the  Douglas  method  of  Hess 
and  Smith.  ' Consider  a differential  element  of  a sink  distribution 
located  at  the  point  P(£,0)  and  having  length  d£,  as  shown  in  the 
sketch.  The  strength  of  the  sink  element  is  qd£.  In  incompressible 
flow,  the  magnitude  of  the  velocity  induced  by  the  differential  sink 
at  a point  P(x,y)  is 


V = qd£/2nr 


(1) 


1/2 


in  which  r is  the  distance  [(x  - + y ] . The  velocity  vector  is 

directed  along  r,  towards  the  sink.  The  horizontal  and  vertical  com- 
ponents of  the  velocity  are  given  by 


u = V cos  0 


qd£  _ (x  - £ ) 
2ttt  r 


(2) 


v = V sin  6 


(3) 
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The  components  of  the  velocity  induced  at  the  point  P(x,y)  by  a 
sink  panel  are  determined  by  integrating  the  contributions  of  all  the 
elements  that  make  up  the  panel.  For  a sink  distribution  of  constant 
strength  per  unit  length  q j , the  horizontal  component  of  velocity  is 
given  by: 


U_H  f 1*- 

271-  J (x  - + 


integration  of  this  expression  yields 


u ■ ln[(x  - <>2  + y2]  |o 


Similarly,  the  vertical  component  of  velocity  is  given  by 


and  integration  yields 


. f 111 

J (x  - O2  + 


3 i -1  , f - xx 

v = 2^  tan  (“T_) 


v = tan 


sx  + y/ 


For  a triangular  sink  distribution,  varying  from  qj  at  x = 0 
to  zero  at  x = ±s,  the  horizontal  and  vertical  components  of  velocity 
are  given  by 
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v=  VI/  7^2  (l  + t)d*  + / ^(1-  *>df 


s 

>/ 


s I 2vr2  s 

o 


integration  of  these  expressions  yields 

q 


u = 


271- 


{s  [tan’1  1 

1 gy 

\ + f -1/ 

-sy  \ 

v x2  + sx  + 

y2 

j i can  i 

x2  - 

sx  + y2  ' 

f + 1)1»( 

(x  + s)2  + 

vV2 

) +<f- 

l)ln 

/(x-  s)2  + 

x2  + y2 

' x2  + y2 

1/2' 


and 


= ^l(f+1)  tan"1(~ ^ 2 

L I S ' X‘  + SX  4- 


+ sx  + y‘ 
1/2 


)+(“  - 1)  tan  1 f — — \ 

' s ' - sv  + / 


- * i„((*  v>2 1 y;)  +u 

s \ x2  + y2  / 


( (x  ~ s> 
V x2  + 


2 + y2 


x2  - sy  + y2  , 
1/2~1  \ 


(8) 


By  suitable  transformation  of  coordinates,  these  expressions  can  be 
used  to  calculate  the  velocity  induced  by  any  panel  at  any  other  point 
in  the  field. 


(7) 
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APPENDIX  II 

SINK  ENTRAINMENT  VELOCITIES 


In  this  Appendix  the  velocities  at  the  midpoint  of  each  sink  panel 
are  derived  in  order  to  relate  the  sink  strength  to  the  entrainment 
velocity.  For  a sink  distribution  having  constant  strength,  q j , over 
the  interval  from  -s  to  s,  the  vertical  component  of  velocity  is  given 
by 


v 


±L  f yd  I 
2?r  J (x  - |)2  + y2 

-S 


(i) 


integration  of  this  expression  gives 


v 


= n 

2* 


2sy 

+ y2  - 


(2) 


Approaching  the  midpoint  of  this  panel,  first  in  the  limit  as  x ->  0, 
then  along  the  centerline  as  y ->  0 yields 

v (0,0)  = (3) 

Thus,  the  sink  strength  per  unit  length  is  given  by  twice  the  entrain- 
ment velocity  qj  = -2Ve. 


For  a triangular  sink  distribution  which  varies  from  qj  at  its  center 
to  zero  at  each  end,  the  vertical  component  of  velocity  is  given  by 
Equation  (8)  of  Appendix  I.  Approaching  the  midpoint  of  the  panel  as 
before,  first  in  the  limit  as  x -*  0,  yields, 

v (0,y)  = ^ ^ ln  ( S g2  7 ) + 2 tan_1 
along  the  centerline.  Then,  as  y -*  0, 

v (0,0)  = (5) 

Thus,  the  strength  of  the  sink  is  given  by  twice  the  entrainment 
velocity,  qj  = -2Ve,  in  this  case  also. 
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